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ABSTRACT 

o 

O ■ We investigate the use of the gravity sensitive neutral sodium (Nal) doublet at 

^ '. 8183 A and 8195 A (Na 8200 A doublet) as an age indicator for M dwarfs. We mea- 

sured the Na doublet equivalent width (EW) in giants, old dwarfs, young dwarfs, and 
candidate members of the /3 Pic moving group using medium resolution spectra. Our 
Na 8200 A doublet EW analysis shows that the feature is useful as an approximate 
I age indicator in M-type dwarfs with ( V-Ks) > 5.0, reliably distinguishing stars older 

C/^ ■ and younger than 100 Myr. A simple derivation of the dependence of the Na EW on 

■ temperature and gravity supports the observational results. An analysis of the effects 

of metallicity show that this youth indicator is best used on samples with similar metal- 
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licity. The age estimation technique presented here becomes useful in a mass regime 



c/3 . where traditional youth indicators are increasingly less reliable, is applicable to other 

^ " alkali lines, and will help identify new- low mass members in other young clusters and 



associations. 



Subject headings: stars: pre-main-sequence, stars: atmospheres, open clusters and as- 
sociations: individual (/3 Pictoris moving group) 

o: 

^SJ . 1. Introduction 

The members of nearby, young moving groups (NYMGs, see Zuckerman &: Song 2004; Torres 
^ . et al. 2008, hereafter T08) are important because: 1) they represent well characterized samples 
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with close, young stars for the study of their astrophysical properties at >3 times the resolution 
possible in the nearest star forming regions, 2) study of NYMG kinematics can infer the recent 
star formation history in the Solar neighborhood, and 3) gas giant planets are still self-luminous 
from gravitational collapse at ages 10 - 100 Myr. This leads to an increased contrast ratio between 
planet and host. Thus, NYMG members are of great interest to planet hunters (see Lafreniere et 
al. 2008; Biller et al. 2010; Wahhaj et al. 2011). However, it is difficult to identify the members of 
spectral types {SpTy) later than M4 because the usual age diagnostics become unreliable. 

West et al. (2008, 2011) showed that the fraction of Ha emitting M dwarfs in the Sloan 
Digital Sky Survey (York et al. 2000) Data Release 7 (Abazajian et al. 2009) spectroscopic sample 
is strongly dependent on SpTy. Nearly all stars later than ~M4 emit Ha, while Ha emission is 
only seen in early- type M dwarfs that are young (<1 Gyr). This is probably a consequence of a 
change in the magnetic dynamo with the onset of fully convective interiors and also applies to more 
energetic emission, such as X-rays, since they are often correlated with Ha (see Covey et al. 2008). 
Li absorption can also serve as an age indicator in low-mass stars, but does not reliably reproduce 
the expected ages of young cluster members when observations arc compared to theory (see Yee 
& Jensen 2010). This discrepancy may be explained by the dependence of Li depletion on stellar 
accretion history (Barraffc & Chabrier 2010). Lithium is also depleted quickly in low-mass stars 
and may be absent by the ages of the NYMGs (Stahler and Palla 2004). 

Gravity sensitive photosphcric lines, such as those of the alkali metals Na and K, can also 
indicate a star's age. It is well known that the neutral sodium (Nal) lines in the optical and near 
infrared are sensitive to photospheric surface gravity in the sense that the stronger the gravity, the 
greater the line equivalent width (EW). As the radius of the star decreases on the approach to 
the main sequence, the photospheric gravity, g = GM^/Rl, increases, and the line widths change 
accordingly. Alkali atoms are easily pressure broadened because of their single valence electron. 
Early studies showed that Nal doublet at --^5900 A {Na D) is useful as a luminosity class indicator 
in late-K and early-M SpTy's (Luytcn 1923: Spinrad 1962). Unfortunately TiO absorption renders 
the Na D EW useless for these purposes beyond SpTy's ~M2. 

Spectroscopic observations of candidate young stars in Taurus, Upper-Scorpius, the Pleiades, 
and the Solar neighborhood reveal that the Nal subordinate doublet at 8183 and 8195 A (Na 
8200 A) is useful for distinguishing young, mid-M and later-type dwarfs from field dwarfs and 
giants (Steele & Jameson 1995, hereafter SJ95; Schiavon et al. 1995; Briceno et al. 1998; Martin 
et al. 2004; Slesnick et al. 2006a, 2006b, 2008; Takagachi & Itoh 2010). These empirical resuhs 
exhibit sufficient EW scatter such that the Nal measurement does not provide a precise estimate 
of age, but does reliably distinguish the stars in young clusters, such as the Pleiades (~125 Myr), 
and billion year old stars in the field (Martin et al. 2010, hereafter MIO). 

Here we concentrate on the doublet at 8200 A (§2) and measure the EW in the spectra of 
old dwarfs, giant standards and candidate members of the P Pic moving group {BMPG, age 10-20 
Myr, T08). We interpret the empirical gravity sensitive behavior of the doublet using calculated. 
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synthetic spectra (§3) and line formation of very optically thick, pressure broadened lines in stellar 
atmospheres (§4). A discussion and summary follow in §5 and §6. 

2. The Na 8200 A Doublet, Sample Selection, and Measurements 

The Na D doublet originates in the 3s ground state of Nal with wavelengths 5890 A (3s^5i/2 — 
3p^^3/2) and 5896 A {3s'^Si/2 - 3p^Pi/2)- The 8200 A doublet {Nal subordinate doublet) is the 
transition between the split 2>p upper level of the Na D lines at 8183 (3p^Pi/2 ~ '^d?^3/2) and 8195 
A {3p'^Pi/2 - 3(fDi/2; Fig. U Bashkin & St oner 1975). Fig. [2] shows a representative sample of 
low-mass dwarf spectra. The Na D doublet becomes difficult to observe in stars with SpTy's later 
than ~M2 because TiO dominates the spectrum near ~5900 A. The 8200 A doublet is flanked by 
TiO and VO absorption bands at shorter and longer wavelengths, and narrow absorption lines of 
photospheric Fe and telluric H2O lie just to the red (e.g. SJ95, their Fig. 2). 

2.1. Samples and Observations 

To gauge the gravity sensitivity of the Na 8200 A doublet we used four samples of stars with 
varying surface gravities. 

High Velocity M Dwarfs: Dwarf stars with large transverse velocities [HV; \tan > 75 km s~^) are 
kinematically consistent with being members of the thick disk population. Thus, they are expected 
to be old (>1 Gyr) and have surface gravity \og{g) ^ 5.0 (Drilling & Landolt 2000, hereafter 
DL2000). We chose 203 solar metahicity, HV K7 - M9 dwarfs from the LSPM catalogs (Lepine & 
Shara 2005a, hereafter LS05a; Lepine & Shara 2005b) for follow up. The metal content of these 
stars was determined using the procedures detailed in Lepine et al. (2007, hereafter L07). These 
stars were observed using several different medium resolution spectrographs on different telescopes 
(S. Lepine, private communication). 

Giants: The surface gravities of red giants lie in the range 1.0 < log{g) < 3.0 (DL2000). We selected 
4 M giants from the Stony Brook/SMARTS Spectral Standards Librar}01§ (Table [1]). The giants 
were observed using the RC Cassegrain Spectrograph on the SMARTS 1.5m telescope at the Cerro 
Tololo Interamerican Observatory (CTIO) by service observers using grating 58 in the first order. 
These settings provide a spectral range 6000 - 9000 A and a resolution of 6.5 A. 

Known f3 Pic Moving Group: Members of the 10-20 Myr old BPMG are expected to have surface 
gravities intermediate to giants and old dwarfs. The radii of these dwarfs are <2 times their main 



^http://www.astro. sunysb.edu/fwalter/SMARTS/spstds. html 

■^SMARTS, the Small and Medium Aperture Telescope System consortium, under contract with the Associated 
Universities for Research in Astronomy (AURA) 
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sequence radius. Therefore, we expect their log(g)'s to be ^0.3 of the field dwarfs. We chose a 
sample of 11 late-type members from T08, Schlieder et al. (2010), and Schlieder et al. (2012, in 
prep.) for follow up (Table[2l). Observations were made during December 2009 using the the Mklll 
spectrograph on the 1.3m McGraw-Hill telescope at the MDM Observatory. We used the 300 line 
mm~^ grating blazed at 8000 A in the first order to produce spectra with a resolution of 6.0 A. 

Young j3 Pic Moving Group Candidates: We studied a subsample of 18 probable young candidates 
[PYCs) selected on the basis of X-ray and UV activity as a youth indicator in Schlieder et al. (2012). 
This subsample is comprised of 2 types of PYCs, those for which the activity based youth indicators 
are reliable {{V-Ks) < 5.0) and those for which activity is an ambiguous youth indicator {{V-Ks) 
> 5.0). We refer to the first subsample as candidates with reliable youth indicators (CWRYs) and 
the latter as candidates with ambiguous youth indicators (CWAYs). These candidates are listed 
in Table [3] and were observed with the known members in 2009 December. 

2.2. Na 8200 A EW Measurements 

All of the MDM spectra were reduced using IRAF. A NeAr arc lamp spectrum was used 
to calibrate the wavelength scale and flux calibration was determined from observations of spec- 
trophotometric standards from Oke (1990). SpTy^s for all dwarf samples were determined using 
the methods described in Lepine et al. (2003). The giant standard spectra were reduced by F. 
Walter using his pipeline procedures (Walter et al. 2004). Giant SpTy^s were taken from the 
Stony Brook/SMARTS Spectral Standards Library. The telluric H2O feature that contaminates 
the 8200 A doublet was corrected in all spectra using custom software written in IDL (S. Lepine, 
priv. comm.). If uncorrected, this feature can increase the measured EW of the doublet by ~10% 
(SJ95) 

The Na 8200 A doublet becomes prominent around SpTy MO (see Fig. [2]). We placed all 
spectra on the same 3 A resolution and measured the doublet in all stars where it was apparent in 
a visual inspection of the spectra. Since the Na 8200 A doublet is unresolved in our spectra, we 
chose a central wavelength (Ac) of 8190 A and integrated the flux in a 22 A region around Ac. To 
normalize the line width, we calculated the average flux in pseudo-continuum regions blueward and 
redward of the doublet. For all dwarf samples these regions were 8149 A - 8169 A and 8236 A - 
8258 A. Slightly different continuum regions were chosen for the giants but Ac remained the same. 
We maintained the same measurement procedure, identical Ac's, and nearly identical continuum 
regions despite the diverse range of instruments used. This minimizes the systematic uncertainty 
introduced by the use of different instruments. Our code allowed widths to be measured on a star 
by star basis while monitoring signal to noise, thereby avoiding particularly noisy measurements. 

Fig. [3] shows the Na 8200 A EW as a function of (V-Kg) for the observed samples. The V 
magnitudes of the dwarf samples were calculated from USNOl.O-B (Monet et al. 2004) B, R, 
and / values following the methods described in LS05a, Kg magnitudes were drawn from 2MASS 
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(Skrutskie et al. 2006). Fig. [3]shows a second order polynomial fit to the HV M dwarf data (dashed 
gray line) with measurement uncertainties. The individual data for the HV M dwarfs is shown in 
Fig. [8] and discussed in §5. The isochrones shown in Fig. [3] are described in §4. For stars having 
( V-Ks) > 5.0 we find a clear distinction between the Na 8200 A line EW of young and old samples. 
A Kolmogorov-Smirnov (KS) test yields a negligibly small probability that the two samples (the 
young dwarfs shown in Fig. [3] and the old M-dwarfs with solar metallicity shown in Fig. [8]) were 
drawn from a common population and differ only by the noise of sampling, of order 10^^. The 2 
known BPMG members with (V-Kg) > 5.0 and several of our CWAYs have EWs which are less 
than those of the old HV dwarfs (see Tables ED . Conversely, a few CWAYs have EWs inconsistent 
with known young stars. This indicates that the Na 8200 A EW is useful in differentiating old, 
active stars from those that are truly young. Johnson and Apps (2009) suggest metallicity has an 
effect on both the EW and the ( V-Kg) color. We discuss the effects in the following section and in 
§5. 



3. Comparison to Models 

We compared our observations to model predictions using Siess et al. (2000, hereafter SDF2000 ji 
and Baraffe et al. (1998, hereafter BCAH98jl| solar metallicity, pre-main sequence (PMS) evolu- 
tionary models combined with synthetic spectra generated using the PHOENIX stellar atmospheres 
code (Hauschildt et al 1999; Rice et al. 2010; hereafter RIO). We obtained effective temperature 
(Tgff) and log{g) data for grids of PMS models with ages 1 Myr - 1 Gyr and masses 0.1 Mq - 0.5 
Mq. To measure the EW of the 8200 A doublet for the models we produced theoretical spectra for 
our {T(.ff, log{g)) grid using PHOENIX. We reduced the spectra to 3 A resolution and measured 
the doublet EW in the same way as the observed dwarf spectra. Fig. H] shows T(,ff =3300K theo- 
retical spectra with 4 log(g)'s which approximately represents the contraction of an 0.2 Mq star to 
the main-sequence. The increase in Na line width with log{g) is apparent. 

In Fig. Owe show trimmed, convolved PHOENIX synthetic spectra compared to observations 
across the entire observed wavelength range. In the top of the figure, observed spectra (black) of old, 
weakly active dwarfs are compared to 1 Gyr old model spectra (red). A likely young, mid-M dwarf 
is compared to a 10 Myr old model (green) in the bottom of the figure. The chosen observed spectra 
have temperatures consistent with the models using the T^fj to SpTy conversion of Kenyon and 
Hartmann (1995, hereafter KH95 jf] for the main sequence stars and the conversion scale of Luhman 
et al. (2003) for the young dwarf (see Table H]). The figure shows that for both main sequence and 
young M dwarfs the observed and model spectra match both the continuum and absorption features 



^See |http://www-astro.ulb.ac.be/~siess/index.html| 
*See http://perso.ens-lyon.fr/isabelle.baraffe/ 

^See also the conversion table by Eric Mamajek at http://www.pas.rochester.edu/~emaniajek 
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well across the observed wavelength range. There are deviations in some of the molecular bands 
due to missing opacity in the models, particularly the TiO band at -^5900 A. This arises from the 
incomplete treatment of complex molecular line lists suffered by the models (Chabrier and Baraffe 
1997; BCAH98). The discrepancy in the TiO band increases as mass decreases in the main sequence 
dwarfs. However, the match is better at this feature in the young M5; this is consistent with the 
results of RIO for high resolution J-band spectra. In their analysis of PHOENIX model fits to late- 
M dwarf spectra in the near-IR they find the molecular features are better represented in young, 
low-gravity models. The comparison in the region around the Na 8200 A doublet is consistent. We 
compared the measured EWs of several observed spectra and their synthetic counterparts. The 
EW differences for the spectra shown in Fig. [5] are +0.4, +0.3, -1.0 and 0.0 A. For several other 
comparisons (not shown) the differences were in the range ±0.2 A. The scatter is less than that 
among the measured EWs as shown in Figs. [3] and [8l and less than the differences attributable 
to gravity at (V-Ks) > 5.0. This simple visual comparison shows that the observed spectra and 
models match well not only in the region of interest around the Na 8200 A doublet, but across 
red, optical wavelengths. Our comparison is motivation for a full analysis of optical model fits to 
observed low-mass dwarfs spanning a range of masses, ages, activity levels, and metallicities. This 
analysis is the focus of another paper (Schlieder et al., in prep.). 

The SDF2000 models provide V and K magnitudes from the empirical Tg/ / to color conversion 
of KH95 (see Table . We checked the conversion from the model K mags to Kg magj^l and found 
the average difference to be ~0.04 mags, which is insignificant in our analysis. In Fig. [3] we show 
only SDF2000 isochrones because the BCAH98 models are consistent when their T^ffs are used 
to interpolate magnitudes from the KH95 scale. We do not show the 1 Myr, 100 Myr, and 1 Gyr 
data point for the 0.1 model because the predicted temperatures at these ages are beyond the 
limits of the KH95 conversion scale. 

A star's metalhcity affects both its 8200 A doublet EW and {V-Ks) color. The Na 8200 A dou- 
blet is saturated, hence its EWcc{NalY^'^ , and to a first approximation, Z^/"^, where Nal is the 
number density of neutral sodium, and Z the metalhcity. We used SDF2000 models to obtain 
colors at Z = 0.01, 0.02, and 0.04 for 0.5 and 0.2 Mq stars at 1 Gyr age. The shifts in the {{V- 
Ks), EW) plane are indicated by the orange vectors in Fig. [3l the dots indicate the extremes of 
metalhcity, twice solar and half solar. The increment in EW is given by its {NalY^^ dependence. 
The increment in [V-Ks) is determined by the change in T(,ff with metalhcity in SDF2000 models. 
Even small differences in metalhcity can produce substantial shifts in the {{V-Kg)-, EW) plane; we 
discuss the effects of metalhcity further in §5. 



®KH95 use Bessel and Brett (1988) near-IR colors for K - M dwarfs, we used the conversion from Kbb to Kg 
given at |http: / /www. cistro.caltech.edu/~jmc / 2mass / v3 /tra nsformations / , See also the magnitude conversion by Eric 
Mamajek at http:/ /www. pas. rochester.edu/~emamajek/memo_BB88_Johnson2MASS_VK_K.txt 
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4. Temperature and Gravity Dependence of Na 8200 A 

In this section we address the fohowing questions: 1) Why does the Na 8200 A doublet become 
prominent in stars later than MO SpTy? 2) How does the strength of the doublet depend on 
photospheric gravity? 3) Why does the EW oi the doublet become useful as a youth indicator for 
{V-Ks)>5.0? 

We used the PHOENIX code to calculate line profiles and evaluate the temperature and gravity 
dependence of the Na doublets in low-mass stars. To calculate the spectra, PHOENIX derives the 
structures of the atmospheres in terms of temperature (T), gas pressure (Pg), electron pressure 
(Pe), and optical depth (r). We calculated two sets of models: 1) Stars of mass 0.7, 0.5, 0.4, 0.2, 
and 0.1 M0 at 1 Gyr age, 2) stars of mass 0.7, 0.2, and 0.1 M0 at ages 1, 10, 50, and 100 Myr, and 
thus model their photospheres as they contract. 

The fraction of Na that remains neutral in the atmospheres of these stars is given by the Saha 
equation 

where Nail, Nal, and are the densities of singly ionized and neutral sodium, and electrons. 
With first ionization potential = 5.1 eV, Pg = NgkT and the ratio of partition functions Zu/Zi = 
1.0/2.0 (Cayrel &; Jugaku 1963), 

Nail T2 59600 , , 

— = 0^33-.-^. (2) 

The second ionization potential of Na is 47 eV; hence at the temperatures of these stars, Nal 1 1 /Nail 
is negligibly small and, to an excellent approximation, the total abundance of Na is Natot = 
Nal + Nail. Fig. [6] shows the results as a function of continuum optical depth in 1 Gyr old stars; 
as Tf,ff of the star decreases an increasingly larger fraction of Na remains neutral. This explains 
the increasing strength of the 8200 A doublet in stars cooler than MO. The increasing fraction of 
neutral Na compensates for the precipitous decrease of the population in the 3p to 3s levels. The 
Boltzmann equation, 

24400 



shows that, 



Nal{3s) gss ' 



Nal{3p) _ 1 1 
Nal{3s) ~ 202 1500' 
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for T = 4000 to 3000 K. The Na D lines would also appear very strong at spectral types cooler 
than MO were it not for their obscuration by TiO, since essentially all Nal is in the ground state, 

iVa/(3s) « ^^Natot oc Pgas- (4) 

To calculate the strength of the 8200 A doublet in absorption we follow the approach of 
Gray (1992, hereafter G92) and write for the increment of opacity in the line, Tune, 

dTline = {f^cont + feline) dx = | 1 H '■ J Hicontdx, (5) 

where is the absorption coefficient in the 8200 A doublet and Kcont is the absorption coefficient 
in the continuum near the lines. The opacity in the line plus continuum is necessarily greater than 
in the continuum alone and the line forms higher in the stellar photosphere than the continuum. 
Since the temperature is lower at the higher levels the line appears in absorption. In low-mass 
stars such as these, Kcont is bound- free absorption of H~. It dominates except for wavelength 
regions where absorption bands such as those of TiO and VO are strong. Absorption by H2O 
becomes important below T^jf ~ 2500K (Auman 1969). The bound- free absorption coefficient of 
if- evaluated at A = 8200 A is (G92, eqns 8.11 and 8.12) 

w(i/-,8200 1) = 0.194^10=^80°/^:^ cx PePgasHT), (6) 

where h{T) contains the temperature dependence. 

The line absorption coefficient is (G92, eqn 11.12 et seq) 

Rine = o^lfl,u4>{l' - ^o), (7) 

where Ni is the population of the lower level, here N{3p), fi^u is the lower to upper level oscillator 
strength, and (j) is the line shape function, normalized such that J (/>(Az^)(i(Az^) = 1, where Ai/ = 
1/ — I'Q and vq is the frequency of the center line. 

The Na D and 8200 A doublets are optically thick (see Fig. H]). The very large line widths, 
several decades larger than the thermal widths, are produced by van der Waals broadening (Curtis 
& Jefferies 1968; Burrows & Volobuyev 2003). It is conventional to designate the van der Waals 
"damping parameter" as (e.g. G92, Eqn. 11.34 et seq.). Far in the wings of the line the first 
term in the Voigt profile suffices, 

0.564 f I V 
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The van der Waals damping term is 

76 = 1.6 X 10^:^ (9) 

Tio 

(Unsold 1955). 

We combine Eqns. (3-4), (6-8), and (9) and find 

Kline OC P^gasf{T), (10) 

where f{T) represents the temperature dependence. Hence, 

^oc%^{|^. (11) 

We used the BCAH98 models for 0.1, 0.2, and 0.7 M© stars at 1, 10, 50, and 100 Myr age 
to determine the dependence of Pgas/Pe on photospheric gravity g at Tconti^-'^t^'m') = 1. Fig. [7] 
shows that Pgas/Pe increases with g in the photosphere and hence with increasing age and explains 
the gravity and temperature dependent nature of the doublet strength. Further, the dependence of 

Pgas/Pe 

on g strengthens as mass decreases, a result consistent with the empirical demonstration 
in Fig. [3l 



5. Discussion 

5.1. Effect of Metallicity on EW and Color 

Fig. [8] shows the Na 8200 A EW as a function of (V-Ks) for 4 samples of M dwarfs with 
increasingly smaller metallicity. In the figure, the solar metallicity, HV M dwarfs are represented 
by black dots. The second order polynomial fit to their distribution (gray, dashed line) is described 
in the caption. The three other samples are K7 - M8 sequences of low metallicity, subdwarf, 
extreme subdwarf, and ultrasubdwarf standards from L07. L07 determined their metallicities from 
the observed TiO/CaH ratio. The ultrasubdwarfs represent the most metal poor M dwarfs in the 
galactic halo. 

The change in EW and color is apparent, lower metallicity leads to smaller Na EWs and 
bluer ( V-Kg) colors, consistent with the predictions in §3. The extreme and ultrasubdwarf samples 
occupy the same area in the Fig. [8] and can be regarded as one sample within measurement errors. 
The KS test applied to the M dwarfs and subdwarfs indicates a high probability, P=0.64, that they 
are drawn from the same population. The test also shows that the extreme and ultrasubdwarfs 
have negligible probabilities, P ~10^^ and 10~^, of being drawn from the same sample as the M 
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dwarfs and subdwarfs respectively. Thus, the application of the Na 8200 A EW to estimate age is 
only reliable for samples of homogeneous metallicity. 

Of the stars observed for the Na 8200 A age indicator analysis, only the HV M dwarfs 
have measured metallicity. However, all of the /3 Pic moving group members and candidates have 
kinematics consistent with the solar neighborhood and most are known to be young. Thus, we 
expect their metallicities to be approximately solar or super-solar. If these dwarfs are very metal 
rich, based on predictions, we expect their positions in Fig. [3] to be systematically shifted to the 
upper right. This is not consistent with the observations. We therefore conclude that the dwarfs 
used in the age indicator analysis represent a homogeneous, approximately solar metallicity sample 
in which stars with systematically small Na 8200 A EWs are young. 

5.2. Application of Results to Other Alkali Lines 

The gravity and temperature sensitivity of the optical Na doublets discussed here is ap- 
plicable to the other Nal lines in Fig. [1] and other alkali atoms and ions with similar electronic 
structure. For example, Greene and Lada (1996) show that the combined EWs of Nal at ~2.2 fim 
and Cal triplet around 2.3 /xm in young stellar objects in p Oph are intermediate to main sequence 
dwarfs and giants, consistent with our results for the Na 8200 A doublet. In addition, McGovern 
et al. (2004) show that KI, Rbl, and Csl in the optical and KI in the near-IR J-band exhibit the 
same gravity sensitivity discussed here. These features are useful as youth diagnostics in the very 
cool photospheres of brown dwarfs. 

6. Summary 

We present an analysis of the Na 8200 A doublet in several stellar samples to show its 
usefulness as an age indicator for dwarfs of mid-M SpTy and later. The samples were chosen to 
be representative of different stages of stellar evolution, from the PMS to evolved giants, which in 
turn trace the evolution of photospheric gravity. We measured the EWs of the 8200 A doublet 
and compared the empirical results to isochrones generated using BCAH98 and SDF2000 evolution 
models and the PHOENIX model atmosphere code. We find that the Na 8200 A doublet is useful 
as an approximate discriminant of youth for stars with (V-Kg) > 5.0, or ~M4 SpTy and later. The 
empirical results are consistent with previous studies and with the theoretically derived dependence 
of the doublet EW on temperature and gravity. The EW of the doublet is useful for mid-M and 
later dwarfs because: 1) it becomes prominent in the coolest dwarfs because of the increased 
abundance of neutral Na and 2) the dependence on gravity is stronger in 0.1 Mq dwarfs than in 
0.7 Mq dwarfs. Analysis of the effects of metallicity show that this youth indicator is best used on 
samples with similar metallicity. The EW technique as we have described it cannot measure an 
age. It can nonetheless discriminate between old M dwarfs in the field and relatively young ones. 



As such, we can expect the technique can be profitably apphed to identify new, low-mass members 
of young clusters and associations. 
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Table 1: Giant Standards 



Name 


a" 


(5« 


V 


Ks 


SpTy 


EW 




(2000.0) 


(2000.0) 


(mag) 


(mag) 




(A) 


HD 1879 


5.768131 


-15.942617 


6.5 


2.3 


M2III 


1.8 


HD 27598 


65.172214 


-16.830004 


7.1 


1.8 


M4III 


1.8 


HD 198026 


311.934330 


-5.0276032 


4.5 


-0.3 


M3III 


1.3 


HD 207076 


326.632640 


-2.212806 


6.8 


-1.7 


M7III 


0.6 



"ICRS epoch J2000.0 RA and Dec in decimal degrees 
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Table 2: BPMG Members 



Name 


a"" 




V 


Ks 


SpTy 


EW 




(2000.0) 


(2000.0) 


(mag) 


(mag) 




(A) 


TYC 1208 468 l'' 


24.414167 


18.592500 


10.7 


6.7 


K7V 


1.6 


HIP 11152'' 


35.860833 


22.735278 


11.6 


7.3 


MIV 


2.1 


AG Tri'' 


36.872083 


30.973611 


10.2 


7.1 


K7V 


1.7 


BD+05 378'^ 


40.357917 


5.988333 


10.2 


7.1 


K7V 


1.7 


PM 104439+3723^ 


70.987083 


37.384250 


13.4 


8.8 


M2V 


3.1 


V 1005 Ori'^ 


74.895000 


1.783611 


10.3 


6.3 


MOV 


2.0 


TYC 1281 1672 1^ 


75.205375 


15.450194 


11.0 


7.6 


K7V 


1.8 


PYC J05019+0108'= 


75.485625 


1.145250 


13.2 


7.7 


M5V 


3.3 


V 1311 Ori'* 


83.018750 


-3.091389 


11.4 


7.0 


M2V 


2.7 


PM 107295+3556^ 


112.379500 


35.933389 


12.1 


7.8 


MIV 


2.5 


PYC J21376+0137'= 


324.417500 


1.620556 


13.6 


7.9 


M5V 


3.6 


BD-13 6424*^ 


353.128750 


-12.264444 


10.9 


6.6 


MOV 


2.2 



"ICRS epoch J2000.0 RA and Dec in decimal degrees 
''Schlieder et al. (2010) 
=Schlieder et al. (2012, in prep.) 
''Torres et al. (2008) 



-17- 



Table 3: Probable Young Candidates'^ 



CWRYs 



Name 


a' 


6' 


V 


Ks 


SpTy 


EW 




(2000.0) 


(2000.0) 


(mag) 


(mag) 




(A) 


PYC J03108+1838 


47.723206 


18.644018 


12.4 


8.8 


MO 


1.8 


TYC 3339 1107 2 


59.333015 


50.855167 


10.6 


7.3 


K7 


1.3 


PYC J03575+2445 


59.391407 


24.752939 


13.0 


8.7 


M2 


2.3 


HIP 47133 


144.066330 


37.529320 


11.1 


7.2 


M2 


2.4 


PYC J10571+0544 


164.297470 


5.748539 


13.2 


9.0 


M2 


2.6 


PYC J11519+0731 


177.986725 


7.523960 


12.7 


7.9 


M3 


4.3 


CWAYs 


PYC J02179+1225'= 


34.483406 


12.424052 


14.4 


9.1 


M4 


2.6 


PYC J03047+2203'= 


46.183666 


22.055891 


15.9 


9.7 


M5 


3.4 


PYC J04116+2504'= 


62.901592 


25.078272 


14.8 


9.4 


M4 


2.5 


PYC J04147+3816'= 


63.691860 


38.267742 


15.7 


9.8 


M5 


3.3 


PYC J07266+1850 


111.673080 


18.842928 


14.2 


9.1 


M4 


3.5 


PYC J08227+0757 


125.697861 


7.954772 


14.6 


9.2 


M4 


4.1 



"Based on strong X-ray/UV activity and (V-Kb) color, see Schlieder et al. (2012) 
''ICRS epoch J2000.0 RA and Dec in decimal degrees 
'^Candidates consistent with youth in Fig. [3] 
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Table 4: SDF2000 Z = 0.02 Model Data" 



Pre-Main Sequence Stars 



0.5 Mq 



Age 




log(3) 


EW 


( V-K)» 


(Myr) 


(K) 




(A) 


(mag) 


1 


3771 


3.5 


1.6 


3.8 


5 


3760 


4.0 


1.7 


3.8 


10 


3764 


4.2 


1.7 


3.8 


50 


3796 


4.5 


1.8 


3.7 


100 


3828 


4.7 


1.9 


3.6 


0.3 Me 


Age 




log(ff) 


EW 


( V-K)" 


(Myr) 


(K) 




(A) 


(mag) 


1 


3360 


3.5 


1.9 


5.3 


5 


3429 


4.0 


1.9 


4.9 


10 


3475 


4.2 


2.0 


4.7 


50 


3531 


4.5 


2.2 


4.4 


100 


3529 


4.7 


2.5 


4.4 


0.2 Mq 


Age 


'^eff 


log(3) 


EW 


( V-K)" 


(Myr) 


(K) 




(A) 


(mag) 


1 


3075 


3.5 


2.3 


7.1 


5 


3235 


4.0 


2.3 


6.2 


10 


3290 


4.2 


2.4 


5.8 


50 


3330 


4.7 


2.9 


5.5 


100 


3317 


4.9 


3.6 


5.6 


0.1 Mq 


Age 




log(3) 


EW 


( V-K)'' 


(Myr) 


(K) 




(A) 


(mag) 


1 


2928 


3.5 






5 


3023 


4.0 


3.0 


7.5 


10 


3075 


4.2 


3.3 


7.1 


50 


3055 


4.7 


5.0 


7.3 


100 


3000 


4.9 






1 Gyr Stars 


SpTyO 




log(g) 


EW 


( V-K)" 




(K) 




(A) 


(mag) 


MO 


3804 


4.9 


1.9 


3.7 


M2 


3669 


5.0 


2.8 


4.5 


M4 


3290 


5.1 


4.2 


5.6 


M6 


2780 


5.3 







"From Siess et al. (2000) models 



From KH95, K is consistent with Kg 
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Fig. 1. — Simplified Nal atomic energy level diagram adapted from Bashkin & Stoner (1975). 
Transition wavelengths are in A . 
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5000 6000 7000 8000 9000 

Wavelength (A) 



Fig. 2. — 3 A resolution spectra of BPMG candidate stars. The spectra show that: 1) The 
unresolved Na D doublet is detected in all spectra but is dominated by strong TiO molecular 
absorption for stars later than about M2. Na D EW measurements would be unreliable for later 
M's. 2) The 8200 A doublet is strong in M dwarfs; its EW is measurable starting at SpTy MO. 
3) Ha emission is common among mid-M dwarfs making it an unreliable youth indicator for those 
SpTy^s. A color version of this figure is available in the online journal. 
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Fig. 3.— Na 8200 A EW as a function of ( V-Kg) for the stellar samples described in the text (see 
key). The dashed gray line is a second order polynomial fit to our HV M dwarf sample (see §5 and 
Fig. [8]). We used the SDF2000 PMS evolution models to connect a stellar model of a specific mass 
and age with a model atmosphere of the corresponding {T^ff, log{g)) to plot isochrones and mass 
tracks. The figure shows that: 1) ( V-Kg) > 5.0 BPMG members and some candidates are clearly 
separated from the HV dwarf sample and are consistent with young model isochrones, 2) Hence 
the 8200 A doublet is a useful age diagnostic for stars of SpTy ~M4 and later (to the right of the 
second, thick, red, dotted line), where activity indicators become unreliable, and 3) for {V-Ks) < 
5.0 the EW is not useful because of quick convergence to the main sequence. We also show the 
region where the Na D line is useful as a luminosity class diagnostic (to the left of the first, red, 
dotted line). The predicted change in EW and color with metallicity is shown in orange. A color 
version of this figure is available in the online journal. 
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6x10" 



4xlO'» — 



2x10" — 



8140 



T,„ = 3300 K 

Log(g) = 3.3 
Log(g) = 3.8 
Log(g) = 4.3 
Log(g) = 4.8 




8160 



0200 A doublel 



8180 8200 
Wavelength (A) 



8220 



8240 



Fig. 4. — Theoretical spectra in the region of the Na 8200 A hnes calculated using the PHOENIX 
stellar atmosphere code binned to a resolution of 3 A. We present models for T^ff = 3300 K over 
a range of log(g). The increase in line width with increasing log(g) is evident. A color version of 
this figure is available in the online journal. 
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0' 
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Young M5 
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log(g)=4.2 





5000 



6000 7000 8000 

Wavelength (A) 



9000 



Fig. 5. — Comparison of observed spectra to models across the observed wavelength range. In the 
top of the figure observed (black) main sequence dwarfs are compared to 1 Gyr old PHOENIX 
model spectra with parameters derived from SDF2000 evolutionary models (red). The models and 
observations generally match well but there are known discrepancies in some molecular bands in 
the models (TiO, ~5900 A) due to incomplete line lists. The bottom of the figure shows a probable 
young, mid-M dwarf compared to a 10 Myr old model spectrum (green). The model is consistent 
with the observation and the agreement at the TiO band is improved; consistent with previous 
results (see RIO). Measured Na 8200 A doublet EW differences between observed spectra and 
their synthetic counterparts are less than the differences attributable to gravity at ( V-Kg) > 5.0. 
The spike around 8700 A in the observed spectra is attributable to instrumental effects. A color 
version of this figure is available in the online journal. 
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rat 
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Fig. 6. — Ratios of Nal/Natot (solid) and H~ /Htot (dashed) as a function of r for 1 Gy old stars. 
The lines show the ratios for the K7- M6 SpTy sequence (bottom line to top line) predicted from 
BCAH98 model data and PHOENIX model atmospheres. In the top of the figure, as Tejj of the 
star decreases an increasingly larger fraction of Na remains neutral. Thus the 8200 A doublet 
increases in strength in stars cooler than MO. In the bottom of the figure, the H ratio remains 
roughly constant with Tg// but increases substantially with increasing r. 
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6.0 



5.5 - 



.S> 5.0 



O 



4.5 - 



4.0 



3.0 



T 1 1 1 1 1 1 1 1 1 r 

Pgas/Pe at t(L2 /xm)=LO 



0.1 Mq 
0.2 Mq 

0.7 




^ g(0.09±0.05) ^ 

10 Myr 



50 Myr 



100 Myr 



1 Myr 



3.5 



4.0 

log(g) 



4.5 



5.0 



Fig. 7. — Pgas/Pe as a function of g for our grid of masses and ages. Pgas/Pe increases with 
age and hence the Na doublets increase in strength. The dependence of the pressure ratio on g 
increases as stellar mass decreases. This is consistent with our empirical result that the EW of the 
8200 A doublet becomes a reliable youth indicator for dwarfs with SpTy later than mid-M. 
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Fig. 8.— EW ofNa 8200 A doublet as a function of ( V-Kg) color for iJF M-dwarfs with samples 
of low-metallicity dwarfs from L07 (see key). The gray, dashed line is a second order polynomial 
fit to the HV M-dwarf data. The uncertainty of the EW measurements across the observed {V- 
Ks) range were estimated by binning the measurements in increments of {V-Kg) and taking the 
standard deviation in each bin. The uncertainty increases toward later-type M dwarfs, probably 
a result of increased noise in fainter targets. As metallicity decreases, dwarfs of like SpTy are 
shifted toward smaller EWs and bluer colors; consistent with predictions. The analysis confirms 
that the Na 8200 A doublet is useful as a youth diagnostic in low-mass M dwarfs for samples with 
homogeneous metallicity. A color version of this figure is available in the online journal. 



